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Abstract 



We examine the gravitational collapse of a non-linear sigma model in 
spherical symmetry. There exists a family of continuously self-similar solu- 
tions parameterized by the coupling constant of the theory. These solutions 
are calculated together with the critical exponents for black hole formation 
of these collapse models. We also find that the sequence of solutions exhibits 
a Hopf-type bifurcation as the continuously self-similar solutions become un- 
stable to perturbations away from self-similarity. 
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I. INTRODUCTION 



The last few years have seen a renewed interest in gravitational collapse, particularly 
with regard to what numerical relativity is able to teach us about the general phenomenon. 
Choptuik's initial discovery of criticality and other behavior strikingly similar to that seen 
in statistical mechanical systems has suggested a deep property of the gravitational field 
equations. 

A good deal of recent work has shown the existence of collapse solutions exactly at the 
threshold of the formation of a black hole for a variety of matter fields. These include both 
real and complex scalar fields vacuum gravity 0|, a perfect fluid and an axion- 

dilaton model from low energy string theory [0]. In each of these models, some common 
behavior emerges. For example, the growth of the black hole mass just off threshold is 
described by a power law relation 



where p is any parameter which can be said to characterize the strength of the initial 
conditions, and p* is the threshold value, i.e., the value for the critical solution. The critical 
exponent 7 is universal within a particular class of matter fields. For example, 7 ~ 0.37 
for the real scalar field, 7 ~ 0.36 for perfect fluid collapse, and 7 ~ 0.2641066 for the 
axion-dilaton (axiodil) system [|T6|. The solutions may also exhibit an echoing behavior in 
that the features of the exactly critical solution are repeated on ever decreasing time and 
length scales. This self-similar behavior of the solutions has been found in both discrete and 
continuous versions. In particular, for vacuum gravity, discrete self-similarity and echoing 
are observed, while in fluid collapse, continuous self-similarity with no echoing emerges. In 
scalar field collapse, both types have been shown to be present. 

The main results of this paper unify the discrete vs continuous self-similarity known 
in the above models. Specifically, we examine a particular non-linear sigma model which 
smoothly interpolates between the complex scalar field model P] and the axion-dilaton model 
as the value of a certain dimensionless coupling constant k varies. We find a family of 
continuously self-similar solutions parametrized by k. Using linear perturbation theory, we 
study the stability of these solutions, and find that the sequence of solutions undergoes a 
bifurcation at a particular value, k ~ 0.0754, where the continuously self-similar solutions 
go from being stable to being unstable. The free complex scalar field (k = 0) is found to be 
on the unstable side of this bifurcation, while the axion-dilaton field = 1) is on the stable 
side. This is in agreement with previous results for both of these matter fields. Further, we 
find that for negative values k < —0.28, the self-similar solutions become ever more unstable 
hinting at the possibility of further bifurcations, more complicated dynamics, and perhaps 
even chaotic behavior in the collapse of these particular models. Since we work only in 
perturbation theory, these latter results are highly tentative, but they suggest the existence 
of more exotic behavior than may have previously been observed. For this reason, full scale 
numerical work on these models would undoubtedly be a very enlightening undertaking. 

Prior to our work, Choptuik and Liebling |n],|T2| recently studied an apparently different 
model, namely Brans-Dicke gravity coupled to a free real scalar field, for various values of 
the dimensionless Brans-Dicke coupling constant, —3/2 < ubd < 00. They use a spherical 




(1) 
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collapse code, and their main result is a change of stability at 00 bd ~ 0. After the con- 
tinuously self-similar solution was found in the collapse of an axion-dilaton field 0, they 
realized that it was their more general Brans-Dicke model for a particular value of ujbd- 
In fact, we find that their Brans-Dicke model is equivalent to some range of our nonlinear 
sigma model (00 > k > 0), with ubd = 00 corresponding to the free complex scalar field, 
and ujbd = —11/8 corresponding to the axion-dilaton field. The bifurcation in stability we 
find here in linear perturbation theory then coincides with the change of stability previously 
found by Choptuik and Liebling; in particular, we agree with their result that, for axion- 
dilaton collapse, the continuously self-similar critical solution is stable and appears to be 
the attractor. The range k < is not present in the Brans-Dicke model, however. 

After the research presented here was completed, but before this paper was posted, 
Hamade, Horne, & Stuart |]T6[ reported numerical and perturbation results on axion-dilaton 
collapse in spherical symmetry. Our results in linear perturbation theory agree with theirs 
with regard to real modes and critical exponents. They also find by a numerical collapse 
code that the continuously self-similar critical solution is stable and is the attractor, in 
agreement with the work of Choptuik and Liebling; this is also consistent with our results 
below on the complex modes. 

The outline of this paper is as follows. In Section II, we give general arguments on what 
kinds of self-coupling of a scalar field may show new critical phenomena in gravitational 
collapse; likely candidates are the non-linear sigma models. In Section III, we introduce the 
particular non-linear sigma model to be studied in this paper, and discuss its relationship 
to matter fields which have been studied previously. Section IV introduces the equations 
of motion, derives their form in the presence of a continuous self-similarity, and sketches 
our numerical approach to solving them. Section V discusses the perturbation of the con- 
tinuously self-similar solutions and the question of stability of these solutions. Section VI 
presents our results and conclusions. 



II. CRITICAL BEHAVIOR AND SELF-INTERACTION 

With the important exception of 0, all the work so far on critical phenomena in gravi- 
tational collapse has assumed spherical symmetry. In spherical symmetry, there is no gravi- 
tational collapse without matter, from Birkhoff 's theorem. Therefore one might expect that 
critical behavior would depend importantly on the model of the matter. Indeed, the critical 
phenomenology and exponents differ among matter models such as real scalar field, ideal 
gas, complex scalar field, axiodil .... However, studying a real scalar field (p, Choptuik found 
that inclusion of a nonlinear interaction term V{(j)) in the action, 

Lmatter = ^V"0V„0 - V {(j)) , (2a) 

V{(f)) = /i^0V2 + A0V4 (2b) 

made no difference in the critical solution itself or in its phenomenology. 

We can understand this result as follows. At least in all known cases, the critical solution 
is either "echoing" (discretely self-similar) or continuously self similar (CSS - admitting a 
homothetic Killing vector field). In either case, by dimensional analysis, the solution cannot 
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depend on any dimensionful parameters. Here we use dimensional analysis appropriate to 
classical general relativity, with a unit of length £ in some system of units where Newton's 
gravitational constant G = 1. A scalar field (real or complex) then has dimensions 
while a Lagrangian must have units It follows that the parameters fi and A above have 
dimensions different from zero; in particular, fi is just the inverse compton wavelength of 
the particle. Since these parameters are dimensionful, the critical solution cannot depend 
on them, consistent with the numerical results. Q 

For this reason we turn attention to a different kind of self-coupling, one which multiplies 
the kinetic term instead of adding to it. The general form is 

^G,j(0^)V>^V«</)-' (3) 

where there are now some number N of scalar fields 0^ (/ = 1 . . . A^), and where Gjj is some 
function of the fields, fixed once and for all to specify the model. The nonlinear functions 
Gij take the form of a Riemannian metric on the internal space of the 0^, the target space. 
Such models are called non-linear sigma models (or "harmonic map" models, as discussed 
by Misner JlOl), and much is known about them in high energy physics, not least because 
they often appear in the low energy limit of superstring theory. By dimensional analysis, 
the scalar fields 0^ are of dimension as is the target space metric G/j. Therefore any 
parameters appearing in Gjj may also be taken as dimensionless, and we can expect the 
critical solution to depend on them. 

What is the simplest nonlinear sigma model we can study? If = 1 then the matter 
action can be reduced to that of a free field by a field redefinition; a 1-dimensional Rieman- 
nian space is always fiat. So the simplest nontrivial value is = 2, wherein the two real 
scalar fields can be grouped into a single complex scalar field 0. For the target space metric, 
the simplest cases are the spaces of constant curvature, namely the 2-sphere, fiat 2-space, 
or the 2-hyperboloid, all with homogeneous metrics. This is the model we shall study. 

III. THE MODEL 

We work with a model defined by the following action 

S^J,*.V^{R-j^^^y (4) 

The complex field F{x^) is a scalar coupled to Einstein gravity with n a real dimensionless 
coupling constant, 

— OO < K < CXD. (5) 



^Choptuik has also tried adding a conformal coupling .^ii0^ to the matter Lagrangian. In contrast, 
^ is dimensionless, so that that critical solution should depend on it. This point deserves more 
investigation. 
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The model given by Eq. (^) is a nonlinear sigma model. As mentioned above, the target 
space of the model is a two-dimensional space of constant curvature. The curvature of this 
internal space is proportional to — k so that the space is hyperbolic for k > and a 2-sphere 
for K < 0. For the particular case k = 1, our model becomes the axion-dilaton (axiodil) field 
f coupled to gravity]^ 

(6) 

I — IT 

It turns out that the value k = 1 is not affected by quantum corrections as it can be 
protected by extended supersymmetry. For k = the model (^) reduces (after a further 
trivial rescaling of the field) to the free complex scalar field coupled to gravity. Thus this 
general model smoothly interpolates between the two particular matter models that we 
have already considered. In fact, for < /t < cxd we find that this nonlinear sigma model is 
equivalent to the model of a massless real scalar field coupled to Brans-Dicke theory. Liebling 
has recently examined this theory using a version of Choptuik's adaptive mesh refinement 
algorithm. He finds behavior qualitatively similar to that found by Choptuik for the real 



scalar field The connection between the two theories can be seen in the relationship 



between the Brans-Dicke coupling constant |T3] ujbd and our constant k 
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ujbd = -77 + < /t < oo. (7) 

Z OK, 

This means that the axion-dilaton model (/t = 1) corresponds to ujbd = —11/8, while the 
free complex scalar field (k = 0) corresponds to ujbd = +00.Q For —00 < k, < the model 
(I) appears not to be equivalent to any Brans-Dicke model; in particular Eq. does not 
apply. The model behaves in a smooth way as k passes through zero. 

Returning to the full model, the field equations in covariant form as derived from the 
action in Eq. are, 

= n TE^ {^aFVbF* + VaF^V^F) , (8a) 

—2kF* 

V^V.F = —-VaFV^F. (8b) 

1 — K\r p 

In this form, these equations are invariant under a global U{1) group of transformations for 
a constant A 

F' = e'^F, -00 < A < 00 (9) 
and which leave the metric unchanged. 



^Notation: We use f here for the axiodil field, instead of r as we did in Q, to avoid confusion 
with logarithmic time coordinate r below. 

^As UJBD — ^ —3/2+, we have k +00; however this may be a funny limit. 
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For K > 0, this model also has an extra global symmetry not present in general relativity, 
namely an S'L(2,1R) symmetry that acts on F, but leaves the spacetime metric invariant; 
this is a classical version of the conjectured SL{2, Z) symmetry of heterotic string theory 
called S'-duality [|l5l. For the axiodil, k = 1, this symmetry acts on f as 

_ af + 6 , . 

where {a,b,c,d) G IR with ad — be = 1, while leaving g^iy invariant. The corresponding 
transformation of F for general k > is 

^ 1 a^F + ^ ^^^^ 



^/K|3*^/KF + a*'' 



where G C with |ap - \f3\^ = 1. 

In the case where k = 0, the extra global symmetry consists of translations in the two 
fiat directions of the target space. Finally, for n < 0, the group of motions on the 2-sphere, 
5*0(3), constitutes the extra global symmetry. 



IV. THE CONTINUOUSLY SELF-SIMILAR SOLUTIONS 

We briefly review the process of setting up the equations such that they are compatible 
with a continuous self-similarity. To begin, we work in spherical symmetry so the metric 
can be taken as 

ds^ = {l + u) [-b^df + dr^] + r^dn^ (12) 

where b{t,r) and u{t,r) are the metric functions. This is essentially Choptuik's metric in 
radial gauge with some minor redefinitions. The timelike coordinate t is chosen so that the 
collapse on the axis of spherical symmetry happens at t = and the metric is regular for 
t < 0. 

We are interested in finding collapsing solutions of our model. In particular we ask 
whether, as in the complex scalar, axiodil, and fluid collapse cases, there exist continuously 
self-similar (CSS) solutions to these equations for arbitrary k. That a spacetime admits a 
continuous self-similarity is described covariantly by the existence of a homothetic Killing 
vector field ^ satisfying 

'^^Cfl'afe = Va6 + Vfe^a = '^Qab, (13) 

where C denotes the Lie derivative. A coordinate system better adapted to our assumption 
of self-similarity involves the coordinates z = —r/t and r = In | — t|. In these coordinates, 
the metric takes the form 

ds"^ = e^^ (il+u 



{b'^ - z^)dT^ + 2zdTdz + dz"^] + z^dn"^) , (14) 
and the homothetic Killing vector is then expressed simply in these coordinates as 
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Cda = dr. 



(15) 



In these coordinates, Eqs. (||) can be written as 



zu — u 



u 



z{u + l) 
ziu + l) 



F\zF' -Py + F'*{zF' - F) 



0^ 



u{u + 1) 



b' 



ub 



= F"A - F + 2zF' + F' 



P b 

ziu - 2) + — (m + 2) - z- 

2; 



+F(7 + 1 - u) + — F*(AF'2 + 2zF'F - F^) 
p 



(16a) 
(16b) 
(16c) 
(16d) 

(16e) 



where the overdot here means djdr and the prime denotes djdz and we define the functions 

^ = b^~z^ P = \-k\F\^. (17) 
The boundary conditions we use are that the solution is regular on the time axis z = 



and on the so called similarity horizon A = b 



0. Regularity on the time axis 2; = at 



the center of spherical symmetry allows us to write the boundary conditions for the metric 
functions 6(r, z) and u(r, z) as 



6(r,0) = l u(r,0) = 



(18) 



The hypersurface defined by A = is where the homothetic Killing vector becomes null. As 
this hypersurface is in the Cauchy development of the initial data, we expect everything to 
be perfectly regular there even though this is a singular point of Eqs. ([TB|). 

The existence of the homothetic Killing vector simplifies these equations somewhat. For 
the general collapse problem without self-similarity, the metric coefficients u and b will be 
functions of z and r, but our assumed symmetry restricts these coefficients to be functions 
of z alone. We could also let the field F be invariant under the action of the vector field 
^, but that would then fail to incorporate the U{1) symmetry which the field equations 
also possess. To allow for more interesting dynamics to occur, we let a U{1) transformation 
accompany the scale transformations (translations in r). Infinitesimally, this amounts to 

C^F = CdaF = iuF. (19) 

This allows us to give the form of F under our assumption of self-similarity as 



F{r,z) 



(20) 



^For completeness, we have included the field equations in the (t, r) coordinates in the appendix. 
However, they are not crucial to our current discussion. 
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The continuously self-similar (CSS) fields are now 

F(r,z)=e-Vo(^) (21a) 
b{T,z)=boiz) (21b) 
u{t, z) = uo{z), (21c) 

where uj is an eigenvalue, determined by solving the field equations. The subscript o that 
we have appended denotes unperturbed values in anticipation of our eventually perturbing 
the exactly self-similar solution. 

Our equations are now just Eqs. ( [T6| ) with the r derivatives of u{z) and b{z) vanishing, 
F and F' being replaced by /o and /q, and F and F being replaced by zcj/q and —uj'^fo 
respectively. Note that with ?io = 0, we can eliminate u'q and we are left with an algebraic 
relation for Uq{z). The equations of motion now reduce to0 

6[, = ^ (22a) 

Ao/^' = r, (~2zujz - z{uo - 2) - %o + 2) - ^^\fo\A (22b) 
\ ^ Po / 

-/o + zivil - no)) - -/o (Ao/^^ + ^Vo') (22c) 

where we have defined 

Ao = bl- z'^ (23a) 
Po = l-«:|/or (23b) 

«o = 4(iN^/o-^/ol' + l/oP)+ (23c) 
Po % 

4 if^izufo - zf,r + f*{tujh - zf,)) (23d) 
Po 

and where the prime now denotes d/dz. 

The boundary conditions at z = for the CSS problem now reduce to 

6o(0) = 1, /o = free real constant, /o(0) = 0, (24) 

where we have used our f/(l) phase symmetry to fix /q as real. We define the value of z 
where Aq vanishes as Z2. As mentioned earlier, we demand regularity at Aq{z2) = and 
this leads to the additional boundary conditions 

bo{z2) = Z2 = free real const /o(-22) = free complex const (25) 

with the constant fQ{z2) being determined by Eq. (^) at the similarity horizon. 



^It is worth pointing out that our notation here is more closely aligned with what we use in our 
paper on the axiodil, k = 1, than the notation in our papers @,14| on the complex scalar field, 
K = 0. 
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Now with the equations and boundary conditions, we can numerically integrate these 
equations. Once we reduce our second order ODE to two first order ODEs and include the 
real eigenvalue uj we have five real equations and five real unknowns. We use our standard 
technique of solving this two-point boundary value problem by shooting with an adaptive 
ODE solver from both boundary points to a point zi in the middle. The free boundary 
values are then found using a Newton's solver for the nonlinear matching conditions. [§] 

We then follow the CSS solution as k varies, and we find that a CSS solution exists for 

-0.60<K<+oo; (26) 

for K = 0, 1 the CSS solution is the same one found in previous work. Our computations 
actually only extend io k < 15, but the behavior is smooth and the CSS solutions seem 
likely to extend all the way to k = oo. On the other hand, our calculations of CSS solutions 
appear to terminate somehow at k ^ —0.60. We are unsure what exactly goes wrong there, 
but we tend to believe that our numerical routine fails and it is not the case that the CSS 
solutions cease to exist for smaller k. It is, however, worth noting that Maison found that his 
sequence of CSS gas collapses terminated at a maximal value k^ax ~ 0.88 where k appears 
in the equation of state for an Eulerian fiuid p = kp. The reason in his case was a change in 
the nature of the eigenvalues associated with the singular sonic point. At k^ax, two of the 
eigenvalues degerate. But we have no evidence that a similar thing occurs here. 

As far as we know, there is only one eigenvalue u possible for the CSS solution for a 
given however we have not looked very carefully for others. 

We also mention that although we describe the spacetime only up to the similarity 
horizon, the spacetime can be continued in these coordinates to z = +oo. This corresponds 
to the spacelike hypersurface t = 0. We expect everything to be regular on this hypersurface 
except at the axis of spherical symmetry since it too is in the Cauchy development of the 
initial data. Thus the apparent singularity in our equations at z = +00 is merely a coordinate 
singularity. By changing coordinates, we can continue the spacetime through t = 0. We do 
not construct this extension here, but as it was possible to make this continuation for the 
complex scalar field and the axiodil cases, we fully expect that such a construction should 
be possible PJ^l- 



V. PERTURBATIONS AND STABILITY 

As interesting as the CSS solutions are, they do not tell us everything we would like 
to know about the gravitational collapse. After all, these are the exactly critical solutions 
p = p* and comprise a set of measure zero in the space of initial conditions of the collapse. To 
reach them, the initial conditions must be tuned with exquisite care. In addition, such things 
as the critical exponents of the black hole scaling relation are found only with information 
gained by collapse slightly away from the critical solution. 

For these reasons, we look to perturbation theory for additional understanding of the 
CSS solutions. It too is not the last word, but it can shed some light on questions of stability 
and in particular allow us to calculate the critical exponents of the black hole growth. 



As described in [14], the very construction of a Choptuon involves stabilization - a bal- 



ancing between subcritical dissipation and supercritical black hole formation with the critical 
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exponent 7 measuring the strength of this black hole/dissipation instabihty. More specifi- 
cally, for initial data close to, but not exactly on the critical solution, the critical solution 
serves as an intermediate attractor with near-critical solutions approaching it but eventually 
running away from it to form a black hole or dissipate the field to infinity. However, in ad- 
dition to this particular instability, we would like to know if there are additional instabilities 
which would rather drive the near-critical solutions completely away from the Choptuon to 
another, perhaps very different, attractor. Thus by appealing to perturbation theory, we 
are looking for both the black hole instability {i.e. the critical exponent) and possibly other 
instabilities indicating the existence of other, stronger attractors. 

So, with the continuously self-similar solutions in hand, we now carry out a linear per- 
turbation analysis of the CSS solutions, still in spherical symmetry. We define the perturbed 
fields as 

b{T,z)^bo{z)+eh{T,z) (27a) 
u{t, z) ^ uo{z) + eui{T, z) (27b) 
F(r,z)^e*--(/o(^)+e/i(r,^)) (27c) 

where again, denotes the 0th order critical solution, 1 denotes the 1st order perturbation, 
uj is the (unique) eigenvalue of the unperturbed equations (which depends on the coupling 
constant k), and where e > is an infinitesimal constant, a measure of how far away the 
solution is from the critical solution in the space of initial conditions. Using Choptuik's 
terminology, we consider the supercritical regime for infinitesimal 

eocp-p*. (28) 

We now perturb the Einstein equations through 1st order in e, to obtain a set of linear 
partial differential equations for the perturbed fields 61, Ui, fi, in the independent variables 
r, z. Following the standard approach, we Fourier transform the 1st order fields with respect 
to the ignorable coordinate r = log(— t). 



u^{a,z) = J e'''^Ui{T,z)dT, (29a) 
bi{a,z) = J e'''^bi{T,z)dT, (29b) 
Ma,z)= [e"'^MT,z)dT; (29c) 



throughout, "will denote such a Fourier transform. The transform coordinate a is in general 
complex. The 1st order field equations now become ordinary differential equations (ODE's) 
in z, and under appropriate boundary conditions, become an eigenvalue problem for a. So- 
lutions of the eigenvalue problem are then normal modes of the critical solution. Generally 
speaking, there will be many different normal modes /i, each belonging to a different eigen- 
value a. Eigenvalues in the lower half plane Imo" < belong to unstable (growing) normal 
modes. Eigenvalues in the upper half a plane correspond to quasi-normal (dying) modes of 
the critical solution. The eigenvalue a is related to the critical exponent by 7 = — l/Imcr. 

We now want to integrate our equations numerically so we need to determine the bound- 
ary conditions. It is important to bear in mind that in addition to solving the equation for 
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fi{a,z), we must also solve the analogous equation for fi{—a*,z)*. Thus, we will have two 
second order ODE's which must be reduced to four first order ODE's, we will have a total of 
six complex equations to integrate. For the perturbation problem, the boundary conditions 
at z = are found to be 

6i(0) = 0, Mi(0) = 0, /((a,0) = 0, /f(-cr*,0)=0, (30) 

/i(cr, 0) = free complex constant, fi{cr*, 0) = free complex constant. (31) 

At the similarity horizon, z = Z2, the boundary conditions are as follows. Both 61(2:2) and 
Ui{z2) are free complex constants. Either /i(o", ^2) or f[{a,Z2) is a free complex constant 
with the other describable in terms of the other boundary conditions at Z2. We chose to let 
f[{<7, Z2) to be free and /i(a", Z2) fixed as this facilitated examining the lower half complex a 
plane. The same is true for the values f[{—a*, Z2)* and /i(— cr*, ^2)*- Counting the eigenvalue 
a, we now have seven pieces of complex boundary data to go with the six complex equations 
we need to integrate. Since the perturbation equations are linear, we expect the solutions 
to scale, so the extra piece of data is merely a reflection of the linearity of the equations. 
Solutions will come in families which will be parameterized by a single complex parameter. 
Thus we have an eigenvalue problem which is well posed and which should yield a discrete 
spectrum of eigenvalues a. 

To solve the 1st order problem we used a Runge-Kutta integrator with adaptive stepsize 
as part of a standard two point shooting method |p, shooting from z = and from both 
boundaries and matching in the middle z = zi. For convenience we solved the 0th order 
system, Eqs. (^), and the 1st order system, Eqs. (^) simultaneously with the same steps 
in z. As discussed elsewhere, the similarity horizon Z2 is a demanding place to enforce a 
boundary condition, and a second order Taylor expansion of the regular solution was used 
for this purpose. 

To solve the 1st order system, we collected all the boundary values but a into a com- 
plex 6-vector X = {fi{a,0) Ji{-a\0)* ,bi{z2), 111(^2) Jii^, Z2)J[{-a* , Z2)*). Because the 
equations are linear, the matching conditions at z = zi are likewise linear in the boundary 
values. A solution is found when the values at zi of (61, Ui, /i(cr), /i(— a*)*, f[{(j), f[{—cr*)*) 
upon integrating from z = match with those found by integrating from z = Z2, for some 
boundary values X. We can express this matching condition 

A{a)X = (32) 

where A{a) is a 6 x 6 complex matrix which is a nonlinear function of cr, constructed 
numerically by integrations of the 1st order equations, Eqs. (|X^), for six linearly independent 
choices of boundary values X. The condition on a for a solution is then 

det A{a) = 0. (33) 

Once a value for a was found that satisfies this condition, the corresponding boundary 
values X were found as a zero eigenvector of the matrix A; these come in one (complex) 
parameter families, as observed above. Solution of Eqs. ([A2| ) with boundary values X yields 
the normal mode itself. Now, A{(7) has been carefully constructed so that it is a complex 
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analytic solution of a. This follows from the fact that all equations leading to A contain a 
but not 0"*, together with some standard theorems about ODE's. Moreover, A{a) has no 
singularities in the closed lower half a plane. These properties allow us to use a number of 
ideas from scattering theory to study det A{a). In particular, there is a theorem for counting 
the number Nc of zeros of det A within any closed contour C in the closed lower half a plane: 

AcArg det A = 27: Nc (34) 

where Arg det A is the phase of det A, and A^Arg det A is the total phase wrap (in radians) 
around the closed contour C, a result similar to Levinson's theorem for counting resonances 
in quantum scattering theory. 

Furthermore, a conjugacy relation holds, 

A*{-a*)=A{a), (35) 

which means that A need only be evaluated for Recr > in the lower half plane. 

The nonlinear equation dety4(cr) = was solved by the secant variant of Newton's 
method [^. The equation being complex-analytic, the 1-complex-dimensional realization of 
the method was used, and it performed well. 

Since our field equations possess gauge invariance due to general coordinate invariance 
and global U{1) phase invariance, some unphysical pure gauge modes will appear at 1st order, 
to the extent that the gauge conditions implicit in our boundary conditions Eqs. (^,0) fail 
to be unique. 

A pure gauge mode arises from an infinitesimal phase rotation cf) — > e*'^</) in the 0th order 
critical solution: 

hi{z) = 0, (36a) 
ui{z) = 0, (36b) 
f\{z)=zfo{z). (36c) 



This gives a time independent solution of Eqs. (|A2|) that satisfies the boundary conditions; 
hence it corresponds to an unphysical mode at a = 

Another pure gauge mode results by adding an infinitesimal constant to time t t + e 
at constant r in the 0th order critical solution. This is possible because our coordinate 
conditions, Eqs. ([T8|) normalize t to proper time along the negative time axis (t < 0, z = 0), 
but the zero of time is not specified. Then the solution is perturbed by 

6i(r, z) =^\r = -{z/t)h\z) = e~^zb'{z), (37a) 
ui{t,z) =^\r = -{z/t)u'{z)=e-^zu'{z), (37b) 
/i(r, z) = e--^(^|,= e-^{-iuU{z) + z^z)). (37c) 

This pure gauge mode has time dependence e"^""^ = e~'^ and so has negative imaginary 
a = —i. 

There are also two more gauge modes which appear as a pair on the real axis. These 
come from the addition of an infinitesimal complex constant, c, to our zeroth order solution: 
F ^ F + ec. The perturbed fields are then 
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bi{T,z) =0, (38a) 
Ui{t,z) =0, (38b) 
/i(r,^) = ce— ^ (38c) 

(38d) 

This mode has a time dependence of e"*'^'^ = e"*"^"^ and so has a = u. Of course, since we 
have A*{—a*) = A{a), the value a = —uj will also solve the equation detA = and be the 
fourth gauge mode.Q Thus there are these four gauge modes in the a plane and no others. 
These modes should appear as numerical solutions, but are unphysical. 

VI. RESULTS AND CONCLUSIONS 

On integrating and solving for the eigenvalues, o"(ft), we found some novel behavior. We 
confirmed the existence of the gauge modes thereby checking the consistency of our method. 
We also found the critical exponent 7(/«) over the range of k values for which we found a 
solution. Figure 1 is a graph of this exponent as a function of the coupling constant. As 
can be seen, the critical exponent for the CSS solution depends strongly on the value of the 
coupling constant. 

In addition, we evaluated detyl(cr) around a large rectangular contour in the lower half 
plane and used Eq. ( p4D to count the zeros lying within. This allowed us to determine if 
there were additional modes in the lower half a plane. Our results were as follows. We did 
find many more modes in the complex a plane. These additional modes are initially in the 
upper half plane for large positive k and approach the real cLXlS clS Kj decreases. Once one of 
these modes crosses the axis into the lower half plane we infer that the leading normal mode 
of the CSS solution has a change of stability. This first occurs aX n ^ 0.0754. We thus have 
the following 

0.0754 <K< +00, CSS stable 
-0.60 < K < 0.0754, CSS unstable 

(39) 

This confirms the original discovery by Choptuik and Liebling of a change of stability at 
ojbd ~ 0; from Eq. (0) the value would be ojbd ~ 0.158. |n],[T^ Note that these results 



are in good agreement with earlier work. The CSS solution for the complex scalar field 
[n = 0) was shown to be unstable by a similar analysis |T^ while the CSS solution for the 



axion-dilaton field (k = 1) was recently shown in to be the attractor in gravitational 



collapse and hence agrees with what we have found here, namely that the solution found in 
is stable. An important question is if the CSS solution becomes unstable at k ~ 0.0754, 
what is the attractor for the collapse. Our conjecture, borne out by collapse calculations of 
Choptuik and Liebling, is that the attractor between < k < 0.0754 {i.e. ujbd ^ 0) is the 



^When we did a similar analysis for the complex scalar field, we were insensitive to these modes 
since we worked with the derivatives of ^(r, z) so that these modes vanished identically. 
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more dynamically interesting discretely self-similar (DSS) or echoing solution analogous to 
the echoing solution originally seen by Choptuik in the collapse of a real scalar field. 

Since everything in our model is smooth at k = 0, as we decrease k below zero, we expect 
the relevant attractor for the collapse to continue to be the echoing solution. However, the 
above mentioned unstable mode turns out not to be the only mode to move into the lower 
half plane. We have evidence for more modes going unstable by k ~ —0.28. We have been 
able to construct these perturbation modes and they appear to be legitimate solutions of 
the perturbation equations and not numerical artifacts. The presence of these additional 
modes suggests that the model becomes ever more unstable, that is more nonlinear, as k 
decreases. So what happens in gravitational collapse as k decreases below —0.28? The 
CSS solution will certainly not be the attractor and the existence of additional unstable 
modes may trigger further bifurcations in the echoing solution. Since our calculations are 
limited to perturbation theory, we can not say this with certainty, but our expectation is 
that the echoing solution will become unstable and bifurcate into an even more dynamically 
complicated solution. One way to determine what happens here with greater assurance 
would be to take a numerical solution for the DSS solution and perform a perturbation 
analysis. That this would be feasible is suggested by Gundlach's results in which he calculates 
the echoing solution as an eigenvalue problem resulting from the assumption of discrete self- 



similarity in the collapse of a real scalar field. [|1^] However, a more direct approach would 
be to perform a full scale numerical collapse calculation in order to understand what is going 
on in this regime. 

In this paper, we have combined a few of the previously studied models of gravitational 
collapse into a single model of a self-coupled complex scalar field. The model is paratmeter- 
ized by a single coupling constant n. In Table 1, we give a summary of some of the key values 
of K. As the value of the coupling constant decreases, the continuously self-similar solution 
which we find undergoes a change in stability. For the regime where the CSS solution is 
unstable, we believe that the attractor for gravitational collapse is an echoing and discretely 
self-similar solution. This change in stability which occus near k = 0.0754 looks like a "Hopf 
Bifurcation", as it is known in the dynamics literature As k continues to decrease, we 
find evidence for additional instabilities in the model, suggesting that there exists at least 
another bifurcation of the collapsing solution. From the lore on other dynamical systems, 
this further conjectural bifurcation might lead to a doubly periodic attractor, or might lead 
to full blown dynamical chaos in gravitational collapse. Additional work will be able to 
determine whether that is indeed the case. 
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Table 1. Range of the model, its relation to the Brans-Dicke model, critical 
exponents, and stability. 



AT 1 • 

N onimear 


BD / scalar 






sigma, K, 




7 


btabihty 01 Lbb 


+00 


-3/2 


< 0.14(?) 


Stable? 


10.0 


-1.4875 


0.1469 


Stable 


1 


-11/8 


0.2641 


Stable 


< 0.0754 


> 0.158 


> 0.373 


Becomes Unstable 





+00 


0.3871 


Unstable 


< -0.28 


n/a 


> 0.435 


Becomes Mucho Unstable? 


< -0.60 


n/a 




(Don't know if CSS exists) 
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APPENDIX A: 



In this appendix, we simply list some of the equations that seemed too cumbersome to 
burden the main portion of the paper with. We have the general collapse equations for our 
model (H) in the usual (t, r) coordinates. 



u 
v! 
b' 




r{u + 1) - 
r{u + 1) 



_p*_p' _)_ _p J?' 
0^ 



uiu + r 



uh 

z 

^^(7^ - T^F) - i.r%F" + 2rhF' + r%' F') 
0^ 



F*{bF 
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IF' 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 



p ■ h 
where ' means d/dt and ' means d/dr, and where 

p= 1-k\F\^. 

The Eqs. (|16D when perturbed as given in Eqs. ( pT]) become our original set as well as 
the following Fourier-transformed first order equations 
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zui + iaul — 
z{uo + 1) 



(zf\ - i{uj + a*) A)* + f* {zf[ - i{uj - a)h) 
z{uq + 1) / ui 2pi 



Pi 



Uo + 



til 



z{uo + 1) 



pI 



2&1, 



1 



+ ^ ((z/^ - tu;fo){zf[ - i{u + a*)A)* + {zf, - zufoTizfl - - a)A)) 



z{uo + 1) 

pI 



Ui 



Uo + 



= -(-^0^1 + uibo) 

z 

= A'Ao + A|2^^ 



(a;-cT) + |/o|' 

Po 



+ - 

z 

4k 



;22(mo-2) + 6^K + 2) 



4kAo , 
"I Jo/o 



(a; - CT + i(l - Mo)) H -fo{izfo + a;/o) 

Po 



+ 



+/r^{w+(^^/o+^/o)'} 



Po 



+6i (260/^' + ^/o7^' + f (^/o + i^/^) + + 2)/4 

[ Po Oo ^ J 



(A2a) 

(A2b) 

(A2c) 
(A2d) 

(A2e) 
(A2f) 

(A2g) 
(A2h) 
(A2i) 

(A2j) 
(A2k) 
(A21) 
(A2m) 
(A2n) 
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FIGURES 

FIG. 1. This is a graph of the critical exponent, 7, of the continuously self-similar solution as 
a function of k. 
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